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We present a trajectory-based semiclassical calculation of the Ehrenfest-time dependence of the 
weak localization correction and the universal conductance fluctuations of a ballistic quantum dot 
with ideal point contacts. While the weak localization correction is proportional to exp(— te/td), 
where te and td are the dot's Ehrenfest time and dwell time, respectively, the variance of the 
conductance is found to be independent of te. The latter is in agreement with numerical simulations 
of the quantum kicked rotator [Tworzydlo et at, Phys. Rev. B 69, 165318 (2004) and Jacquod and 
Sukhorukov, Phys. Rev. Lett. 92, 116801 (2004)]. 
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I. INTRODUCTION 

Weak localization and universal conductance fluctu- 
ations are manifestations of quantum interference on 
quantum transport. 1 How these effects appear in a ballis- 
tic quantum dot, coupled to source and drain reservoirs 
via ideal point contacts, is believed to depend on the ra- 
tio of the Ehrenfest time and the mean dwell time td. 2 
The Ehrenfest time te is the time it takes for two clas- 
sical trajectories, initially a Fermi wavelength apart, to 
diverge and be separated by a distance comparable to the 
system size^^ Trajectories need to remain inside the dot 
during at least a time te in order to contribute to weak 
localization or universal conductance fluctuations. 

For the weak localization correction SG, semiclas- 
sical theory predicts that 5G oc exp(— te/td) ) 2 i 5 i 6 i 7 i 8 
which has been verified using accurate numerical sim- 
ulations of the quantum kicked rotator^ — a "quantum 
dot" with "stroboscopic" dynamics^ For universal con- 
ductance fluctuations, numerical simulations found no 
Ehrenfest-time dependence of varGii2*ii This remark- 
able result remained unexplained because of the absence 
of a semiclassical theory of universal conductance fluc- 
tuations in ballistic quantum dots. It is the goal of this 
article to report such a theory and offer a microscopic ex- 
planation for the absence of Ehrenfest-time dependence 
of conductance fluctuations seen in the numerical simu- 
lations. 

Two theoretical approaches have been taken to address 
quantum transport in ballistic quantum dots. One is the 
theory of Aleiner and Larking which considers quantum 
corrections to the ballistic analogue of the "diffuson" and 
"cooperon" propagators that play a central role in the di- 
agrammatic perturbation theory for disordered conduc- 
tors. The other approach is based on an expression re- 
lating the dot's scattering matrix to a sum over classical 
trajectories connecting the two point contacts. 3,2 In all 
cases where both approaches have been used to calculate 
the same observable, the results have been the same. In 
this article, we we use the trajectory-based approach. 

In the trajectory-based semiclassical approach, the 
dot's conductance G is calculated from the Landauer for- 
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where T is the total transmission of the dot. The trans- 
mission T is expressed as a double sum over classical 
trajectories a, (3 that connect the two point contacts^ 
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Here A a and Ap are stability amplitudes and S a and 
Sp are the classical actions of the two trajectories a and 
[3. The classical trajectories a and (3 start with initial 
transverse momentum compatible with the same mode n 
in the left contact and end with transverse momentum 
compatible with the same mode m in the right contact. 
The modes in each contact have quantized transverse mo- 
mentum 



p±(m) — inhm/Wj , m = 1, 



(3) 



where Wj is the width of the contact, Nj the number of 
modes in the contact, and the subscript j — 1,2 refers to 
the left and right contacts, respectively. 

Weak localization is the small negative correction 5G 
to the ensemble average of the dot's conductance G aris- 
ing from quantum interference. (The ensemble average 
is taken with respect to small variations of the shape of 
the quantum dot or the Fermi energy.) Two different tra- 
jectories a and (3 contribute to weak localization if their 
action difference S a — Sp is of order Ti. Pairs of trajec- 
tories a and (3 that contribute to weak localization are 
shown schematically in the top panel of Fig. ^ The two 
trajectories are almost equal, except for a stretch where 
(3 has a small-angle self-intersection and a has a small- 
angle avoided self-intersection. Such pairs of trajectories 
were originally pointed out by Aleiner and Larking they 
were first investigated in the trajectory-based formalism 
by Sieber and Richter>i4 The action difference between 
the two trajectories is of order h if the duration of the self- 
encounter is of order of the Ehrenfest time. The prob- 
ability that the trajectories do not escape through the 
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FIG. 1: Top: schematic picture of a pair of interfering trajec- 
tories a and j3 that give rise to the weak localization correction 
to the conductance. Bottom left and bottom right: schematic 
picture of quadruples of interfering trajectories ai, Pi, a.i, 
and P2 that contribute to conductance fluctuations. The true 
trajectories inside the dot are piecewise straight, with specu- 
lar reflections off the dot's boundary. The small-angle (self) 
encounters are shown thick. In the bottom right panel, the 
trajectories j3\ and ai have one more revolution around a pe- 
riodic trajectory (shown dotted) than their partners a\ and 
/3 2 . 



point contacts during the duration of the encounter is 
exp(— te/td), hence the suppression 5G oc exp(— te/td) 
of the weak localization correction for large Ehrenfest 
times mentioned previouslyi^*Si^ 

Since the transmission T is a double sum over classi- 
cal trajectories a and /3, the variance of the conductance 
is expressed as a quadruple sum over classical trajecto- 
ries a%, fix, «2, and /?2- Quadruples of trajectories for 
which the individual action differences S ai — Sp t and 
S a2 — Sp 2 are large (so that they do not contribute to 
the average conductance), but the total action difference 
S ai + S a2 — Sp ± — Sp 2 is small (of order Ti) contribute 
to the conductance fluctuations. The two distinct classes 
of quadruples of trajectories that contribute to varG are 
shown in the bottom panels of Fig. Both classes have 
their counterpart in the theory of conductance fluctua- 
tions in diffusive conductors 

The trajectories in the bottom left panel of Fig. ^ are 
the ones commonly associated with conductance fluctu- 
ations in ballistic quantum dotsA&iLi&iii (see, however, 
Ref. |2(J for an exception) . They undergo two small-angle 
encounters. The trajectories ax and fix, and ct2 and 02 
are pairwise equal before the first encounter and after the 
last encounter. Between the encounters ax and P2, and 
02 and f3x are pairwise equal. The total action difference 
S ai + S a2 — Sp ± — Sf3 2 is of order h if the duration of each 




trajectory 

FIG. 2: If the phase space distance between the trajectories 1 
and 2 is of order Ti at the point when they arrive at or depart 
from the periodic trajectory, their motion is correlated for a 
time ~ te while away from the periodic trajectory. The con- 
tribution of such trajectories to the conductance fluctuations 
is suppressed if te S> td. After the removal of trajectories 
with correlated motion away from the periodic trajectory the 
net contribution of the trajectories shown in the bottom right 
panel of Fig. to the conductance fluctuations zero is finite. 



encounter > te. Since the survival probability during 
each encounter is exp(— te/td), the contribution of these 
trajectories to varG decreases oc exp(— 2te/td) for large 

TE- 

There is no exponential suppression at large Ehrenfest 
times for the trajectories shown in the bottom right panel 
of Fig. 2] Here, the trajectories (3x and a2 differ from ax 
and 02, respectively, by one extra revolution around the 
same periodic trajectory (dotted). The additional revo- 
lution around a periodic trajectory ensures that the indi- 
vidual action differences S ai —Sp 2 and S a2 —Sp 2 are large, 
whereas the total action difference S a 1 +S a . 2 —S/3 1 —Sp 2 is 
small if the trajectories ax and P2 spend at least a time 
te near the periodic trajectory. For a dot with mean 
dwell time td, the typical period of a periodic trajectory 
(weighed with the square of the stability amplitudes) is 
of order td. This means that the trajectories ax, 0x, <%2, 
and /?2 need to wind many times around the periodic tra- 
jectory if they are to contribute to conductance fluctua- 
tions if te td- However, since the survival probability 
of the trajectories ~ exp(— t p /td) depends on the period 
t p of the periodic trajectory only, and not on the actual 
time ~ te spent inside the quantum dot, the contribu- 
tion to the conductance fluctuations is not suppressed 
exponentially if te 3> td- 

This is not the full story, however. Periodic trajec- 
tories are related to density of states fluctuations via 
Gutzwiller's trace formulajSi so that trajectories of the 
type shown in the bottom right panel of Fig. ^ represent 
the impact of density-of-states fluctuations on the con- 
ductance. However, for chaotic quantum dots with ideal 
point contacts, the density of states and the conductance 
are known to be statistically independent^ at least in 
the regime te <C td in which random matrix theory is 
valid. In other words, there is no Einstein relation for the 
conductance of a chaotic quantum dot with ideal point 
contacts if te <C td . The same conclusion can be drawn 
based on the trajectory-based approach: Indeed, if one 
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sums the contribution of trajectories of the type shown 
in the bottom right panel of Fig. over all initial and 
final positions where the trajectories merge with or de- 
part from the periodic trajectory, the net contribution to 
var G vanishes. 

The key observation that explains the existence of 
Ehrenfest-time independent conductance fluctuations in 
ballistic quantum dots is that the cancellation that is 
responsible for the vanishing of the density-of-states con- 
tribution to the conductance fluctuations is lifted if the 
dot's Ehrenfest time te is large in comparison to the 
mean dwell time ttj. This is because correlations be- 
tween the trajectories shown in the bottom right panel 
of Fig.^persist away from the periodic trajectory if they 
approach the periodic trajectory and/or leave it at very 
close phase space points, as shown in Fig.|3 These addi- 
tional correlations, which last up to a time rg, suppress 
the contribution of these trajectories to the conductance 
fluctuations if te ^> td, even if the two trajectories enter 
or exit through the same contact. After removing the 
trajectories with these additional correlations from the 
summation over all trajectories of the type shown in the 
bottom right panel of Fig. ^ the summation over the 
remaining trajectories no longer vanishes. The persis- 
tence of conductance fluctuations if te 3> td thus can be 
attributed to all trajectories of the type shown in the bot- 
tom right panel of Fig. that do not arrive at or depart 
from the periodic trajectory at close phase space points. 

These qualitative arguments will be supported by the 
semiclassical calculations presented in the next three sec- 
tions. We use the trajectory-based semiclassical formal- 
ism in the formulation developed in a series of works by 
Haake and collaboratorsi 23 ' 24 i 25 i 26 In Ref. |U Heusler et 
al. show how this formalism is applied to the calculation 
of the weak localization correction, correcti ng tw o cancel- 
ing mistakes in the earlier theories of Refs. 1511271 (see also 
Refs. 00)- Although developed for the limit te <C td, 
the calculation of Ref. is readily extended to include 
Ehrenfest-time dependences. The beginning of our calcu- 
lation follows that of Heusler et al., but we take a different 
classical limit at the end of the calculation. Heusler et al. 
take the classical limit Ti — > while keeping the number 
of channels iVi and N2 in the two point contacts fixed. If 
the classical limit is taken this way, the ratio te/td — > 0, 
so that the Ehrenfest-time dependence of the weak local- 
ization correction SG and the conductance variance var G 
are lost. In order to preserve the Ehrenfest-time depen- 
dences of 5G and varG, we take the limit h — ► while 
keeping the ratio te/td fixed. For this classical limit, 
both the channel numbers Ni and N2 and the dwell time 
td diverge, although the divergence of the dwell time is 
only logarithmic in Ti. The divergence of the channel 
numbers is of no concern for a calculation of the weak lo- 
calization correction or conductance fluctuations, since, if 
N\ and N2 are large, both quantities depend on the ratio 
Ni /N2 only. The divergence of the dwell time suppresses 
non-universal contributions to the quantum interference 
corrections. Finally, in this limit, interference of trajec- 




FIG. 3: Detail of the small-angle self-encounter and the def- 
initions of the various times used in the calculations. The 
beginning and end of the encounter region are marked by 
solid thin lines, the location of the Poincare surface of sec- 
tion is marked by a dashed line. Each trajectory passes the 
Poincare surface of section twice. Left: Small-angle self en- 
counter fully inside the dot. Right: Small-angle self-encounter 
that touches the lead opening. One end of the encounter re- 
gion is marked by a solid thin line, the other end is the lead 
opening. 



tories with more than one small-angle self encounter (for 
weak localization) or with more than two small-angle en- 
counters (for conductance fluctuations) can be neglected 
since the contribution of such trajectories to the average 
conductance is of order l/(Ni + N2) or smaller^ 

Before we discuss our calculation of the Ehrenfest-time 
dependence of universal conductance fluctuations, we re- 
view the trajectory-based calculation of the weak local- 
ization correction to the conductance G. This allows us 
to introduce the necessary formalism in a calculation that 
is less complex than the calculation of conductance fluc- 
tuations. 



II. WEAK LOCALIZATION 

Without quantum interference, the ensemble average 
(G) = {2e 2 /h)N 1 N 2 /(N 1 + N 2 ), where Ni and N 2 are 
the numbers of propagating channels in the dot's left and 
right point contacts, respectively. We are interested in 
the small deviation SG of (G) from its classical average. 

Starting point of the calculation is the Landauer for- 
mula flj together with the semiclassical expression 
for the dot's transmission. The pairs of classical trajec- 
tories a and j3 that contribute to weak localization are 
repeated schematically in the left panel of Fig. [31 together 
with the definitions of the various times used in the cal- 
culation below. Following Refill (and referring there for 
details) we find that the weak localization correction to 
the transmission is equal to 

ST = N * N > ( FT f ^Le-^A 

POO PC p isu/h— i e nc/TD 

x / d T3 e- Ta/To / dsdu - . (4) 

JO J-c ZTThtenc 
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Here T\ and t 2 are the durations of the parts of the trajec- 
tories between the encounter and the lead openings and 
t 3 is the duration of the loop, see Fig. The prefactor 
NiN 2 /(Ni + N 2 ) 2 consists of a factor Ni from the free 
sum over the incoming channel in the left contact, a fac- 
tor N2/(Ni+N2) from the probability to escape through 
the right contact, and a factor l/(Ni + N 2 ) from the 
small-angle phase space encounter^ The phase space co- 
ordinates u and s parameterize the distance between the 
two parts of one of the interfering trajectories at a refer- 
ence Poincare surface of section during the self-encounter, 
along unstable and stable directions in phase space, re- 
spectively, see Fig. [3J The product su in the exponent of 
Eq. (0J is the difference of the actions of the trajectories 
a and {3m&2& The integration domain is the set of phase 
space coordinates s and u for which the distance between 
the two trajectories is smaller than a cut-off c, i.e., \s\, 
\u\ < c. The cut-off c is a classical scale chosen small 
enough that the classical dynamics can be linearized on 
phase space distances below c. The precise value of the 
cut-off is not relevant in the limit h — ► 0. The time t enc 
is the duration of the encounter, i.e., the time that the 
two segments of the trajectories a and (3 are within a 
phase-space distance c, 



1 , c 1c 
-In — + - In— , 

A s A \u\ 



(5) 



where A is the Lyapunov exponent of the classical mo- 
tion inside the quantum dot. It appears in the denomi- 
nator in Eq. Q in order to cancel a spurious contribu- 
tion to the integral from the freedom to choose the refer- 
ence point inside the encounter regioni^Si The factors 
exp(— tenc/To) and exp(— t^/t-q) are the probabilities not 
to escape during the encounter and during the loop seg- 
ment, respectively. Although the trajectory traverses the 
encounter twice, the escape probability involves only one 
encounter durationAS 

The calculation of ST closely follows the general prin- 
ciple outlined in appendix D3 of Ref. 0- Limiting the 
integration domain to positive u, we first rewrite Eq. Q 
as 



0} becomes 



ST 
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{N 1 +N 2 ) 2 ir J 
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dxx 1 ' Atd cos(xr) 



{Ni + N 2 y-K 



Atd sin r 



SIM 



(8) 



The term proportional to sinr is a rapidly oscillating 
function of Planck's constant and is discarded in the clas- 
sical limit. Writing the remaining term in terms of the 
Ehrenfest time, 



1 1 c 2 

r E = -lnr = T hi — . 
A A n 



(9) 



we find 



8T = 



NiN 2 



(iVi + N 2 y it 



d^lATD^ (10) 



The limit h — > at a fixed ratio of Ehrenfest time and 



dwell time consists of sending both 



2 /% — * oo and 



Atd — > oo. These limits can be taken independently in 
Eq. HDJ. We then find 



ST 



N X N 2 



-tb/td 



(11) 



The exponential dependence of Eq. is in agreement 
with previous calculations of the Ehrenfest-time depen- 
dence of weak localizationi£i£iL& Note that the appear- 
ance of the classical phase-space cut-off c in the definition 
of the Ehrenfest time does not affect the final result: In 
the limit h — -> at fixed te/tq, the dwell time td — > oo, 
which removes any c-dependence from the final expres- 
sions. 

Instead of calculating the weak localization correction 
to the transmission, one may also calculate the quantum 
interference correction to the ensemble averaged reflec- 
tion. Unitarity relates the total reflection Ri off contact 
i to the total transmission T , 



ST = 



2NiN 2 Ty 

(JVi + JVa) 2 



du I ds 



e *™ c / TD cos(su/H) 



2nht r 



(0) 



We then perform the variable change 



c/a, 



cxa. 



(7) 



With the new integration variables, the integration do- 
main is —1 < x < 1 and 1 < a < l/\x\. Further, 



Ri = Ni-T, * = 1,2. 



(12) 



Although this implies SRi — 8R 2 = —ST, it remains 
instructive to verify this result explicitly from the semi- 
classical formalism. 

The semiclassical formula for the total reflection Ri is 
the same as Eq. for the total transmission T, but 
with a double sum over trajectories that connect contact 
i to itself. The quantum correction to Ri consists of two 
parts: the counterpart 5R-1 of the weak localization cor- 
rection to the transmission T, which involves encounters 
in the interior of the quantum dot, and an extra quan- 

(2) 



A 1 In ( 1 / 1 a; | ) . Hence, with r — c 2 /h, the integral turn correction 8 R\ from encounters that touch the lead 
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opening. The calculation of the first correction to Ri only 
differs from the calculation of ST in the replacement of 
iVijV 2 by Nf , i = 1, 2, so that 



5R 



(i) 



N? 



(iVi + N 2 y 



-t e /t d 



(13) 



(2) 

The calculation of the reflection correction SR\ 
from encounters that touch the lead opening is a little 
different This correction is usually referred to as 
the 'coherent backscattering' correction to reflection. In 
the limit te -Ctd, coherent backscattering can be calcu- 
lated using the 'diagonal approximation' for the double 
sum over trajectories in Eq. ^M^^ In the diagonal 
approximation, only trajectories a and (3 that are identi- 
cal up to time reversal are kept. If te ^ td, the diagonal 
approximation fails, however, and a full summation over 
families of trajectories very similar to the trajectory sums 
for weak localization is called forA^i 

Since the trajectories a and f3 have the same perpen- 
dicular component of the momentum at the lead opening 
upon entrance as well as upon exit, the optimal choice 
of phase space coordinates for an encounter that touches 
the lead opening is a phase space coordinate v that rep- 
resents the perpendicular momentum at the contact, to- 
gether with the unstable phase space direction u (taken 
with respect to motion away from the contact). Note 
that v is a meaningful coordinate for a Poincare surface 
of section in an encounter that touches the lead opening, 
because all trajectories piercing the Poincare surface of 
section exit the quantum dot together. The stable phase 
space coordinate s, which was used for small- angle self- 
encounters that contribute to weak localization, is not a 
good choice here, because in general the trajectories a 
and j3 have different s if the self encounter touches the 
lead opening. Following Refs. 12^12^ . we normalize v such 
that the cross section volume element in phase space is 
dudv. With this normalization, the phase space coordi- 
nates (u, v) are uniformly distributed for ergodic motion. 
As in the case of the calculation of the weak localization 
correction ST, we use the first passage of the trajectory 
a through the encounter region as a reference, consider 
a Poincare surface of section at an arbitrary point dur- 
ing the encounter, and denote the phase space coordi- 
nates at which a cuts through the Poincare surface of 
section a second time by (u,v). The action difference 
AS is the phase space area enclosed by the four segments 
of the trajectories a and (3 involved in the interference 
correction s^ 8 ' 2 - Recalling that a and its partner trajec- 
tory (3 have perpendicular momenta compatible with the 
same modes in the lead opening, the phase space coordi- 
nates of the first and second passages of (3 are (u, 0) and 
(0,v), respectively, hence AS* = uv. Note that the en- 
closed phase space areas are conserved along the motion 
of the trajectories. In particular, this means that the ac- 
tion difference AS" is independent of where the Poincare 
surface of section is chosen and that the coordinate v 
scales oc exp(— Xt) upon moving away from the contacts. 



With this, we find that the coherent backscattering 
correction reflection becomes 



SR 



(2) 



Ni 



dT*e- T3/TD 



{N x + N 2 ) 

"A -1 ln(c/|u|) j T giUTt/fi-teno/TD 

td 2nht cnc 



dvdu 

, (14) 



where r is the time needed to go between the Poincare 
surface of section and the lead opening, see Fig. [3] The 
prefactor Ni arises from the summation over the incom- 
ing transverse channel. Note that there is no additional 
factor Ni/(N± + N 2 ) [as in the case of an encounter that 
resides in the interior of the quantum dot, cf. Eq. Q], be- 
cause the probability of escape through contact i is unity 
for encounters that touch the lead opening, i = 1,2. In 
Eq. p4fl the encounter time is a function of r and u, 



i e nc(T, u) = t + X 1 ln(c/|w|). 



(15) 



The condition r < A -1 ln(c/|u|) in Eq. (|14fl ensures that 
the encounter touches the lead opening. (The precise 
value of the cut-off c is unimportant in the limit fi — > 0, 
as before.) In order to perform the integrations in Eq. 
iJHJl, we take u to be positive and perform the variable 
change 



r + A 1 ln(c/|w|), u = c/a, v = cxa. 



(16) 



With the new integration variables, the integration do- 
main is— 1 < at < 1, 1<<7< e Ar , and < t' < 
A -1 ln(l/|a:|). Further, t enc = r'. Hence, with r = c 2 /Ti, 
the integral l|14(l becomes 



SR. 



(2) 



dx 
1 Jo 



\- x U{x/\x\) dT , 



-t' /td 



tt(JVi +N 2 ) 

x / — - cos(rx) 
Ji °r 



(17) 



Neglecting terms proportional to sinr, the remaining in- 
tegrals are the same as for the calculation of ST, and one 
finds the result 7 ! 30 ? 31 



SR) 



(2) _ 



Ni 



m+N 2 



^-t e /t u 



Adding Eqs. {T3J) and (JTSJl, one finds SR 
required by unitarity. 



(18) 
-ST, as is 



III. CONDUCTANCE FLUCTUATIONS 

The same theoretical framework can be used to calcu- 
late universal conductance fluctuations. Technically, it 
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FIG. 4: Two configurations of interfering trajectories that 
contribute to the covariance cov(i?i, R2) of reflection coeffi- 
cients for the left and right contacts. The encounter regions, 
segments of the trajectories for which the phase space distance 
between the trajectories 1 and 2 is less than a classical cut-off 
c, are shown thick. In the top configuration, the trajectories 
1 and 2 have two consecutive well-separated encounters. In 
the bottom panel, the trajectories Q2 and /3i have one more 
revolution around a periodic trajectory (shown dotted) than 
their counterparts ai and P2. 




FIG. 5: Two examples of interfering trajectories with overlap- 
ping encounters. The encounter regions are shown thick. In 
both panels, the trajectories /3i and «2 have one more revolu- 
tion around a closed periodic trajectory (shown dotted) than 
their partners ot\ and P2. The top panel shows two encounters 
that overlap at one end. This is the "three encounter" of Ref. 
|2.1 . The bottom panel shows two encounters overlapping at 
both ends, so that the encounter region forms a closed loop. 
Although encounters that themselves form a closed loop do 
not contribute to the conductance variance if te < td, their 
contribution is essential if te ~S> td. 



is most convenient to calculate the covariance of reflec- 
tion coefficients i?i and R2 for reflection from the left 
and right point contacts, because it avoids the necessity 
of dealing with encounters that touch the two point con- 
tacts. (The case of encounters that touch the lead open- 
ings will be discussed at the end of this section.) The 
covariance of the reflection coefficients is directly related 
to the conductance variance, 

v&rG=(jp\ va,vT=(^p\ cov (i? 1; R 2 ). (19) 

The reflections R\ and R2 arc expressed in terms of a 
double sum over classical trajectories at and ft, i = 1,2, 
connecting each contact to itself, similar to the trajectory 
sum of Eq. J2J) for the dot's transmission. The covariance 
cov (Ri, R2) then becomes a quadruple sum over classical 
trajectories a\, ft, and [32- There are two distinct 
configurations of four interfering trajectories that con- 
tribute to cov (i?i, R 2 ), see Fig. 0] These are the coun- 
terparts for reflection of the trajectories shown in the two 
bottom panels of Fig. ^ We refer to them as "type a" and 
"type b" interfering trajectories^ Other possible config- 
urations of interfering trajectories will give contributions 



to var G smaller by a power of Ni + N2 and need not be 
considered in the limit h — > at fixed te/td we consider 
here^ The configuration of Fig.^J) is not left-right sym- 
metric and acquires an extra factor two. Denoting the 
contributions from interfering trajectories of types a and 
b as A and B, respectively, we then have 

/2e 2 \ 2 

varG = 2 ( — J (A + 2B). (20) 

The prefactor 2 in Eq. I|20l) only appears in the pres- 
ence of time-reversal symmetry. It accounts for the con- 
figurations of interfering trajectories obtained by time- 
reversing the trajectories originating from the left con- 
tact but not the trajectories originating from the right 
contact. 

The primary distinction between the types a and b is 
that the trajectories of the former type do not have seg- 
ments that are close to a periodic trajectory, whereas in 
the latter case they do. (The periodic trajectory is the 
dotted loop in the right panel of Fig. 0] Here "close" 
means that a segment of the trajectories a or /? can be 
deformed into a periodic trajectory within the region of 
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phase space in which the chaotic dynamics in the quan- 
tum dot can be linearized.) This means that B contains 
all conductance fluctuations that are tied to density of 
states fluctuations, whereas A represents the density-of- 
states independent fluctuations of the conductance. In 
the limit te <C td the conductance and the density of 
states are statistically independent^ so that we expect 
that B = in that limit. 

Although the encounters as drawn in Fig. 0] do not 
overlap, they can overlap in principle. In fact, such 
overlaps are essential for a theory of universal conduc- 
tance fluctuations if the Ehrenfest time is larger than 
the dwell time and, hence, comparable to the total path 
length. Two examples of overlapping encounters are 
shown in Fig.0 The top panel of Fig.0shows the "three- 
encounter" of Ref. |23j The bottom panel shows a more 
complicated configuration in which the two encounters 
overlap both at their beginning and at their end, so that 
the encounter region forms a closed loop. In principle, 
overlapping encounters can arise from bringing the two 
encounters of Fig. @Ji close together along one of the two 
solid trajectories. However, as soon as the encounters of 
Fig-UK overlap, one of the trajectories involved has a seg- 
ment close to a periodic trajectory. Therefore, according 
to our definitions of the covariance contributions A and 
B, the two encounters in the configuration of type a do 
not overlap. 

Calculation of the contribution A of interfering tra- 
jectories of type a is straightforward since the quantum 
interference correction from non-overlapping encounters 
factorizes^ 

A= (W'l'-**/* (21) 

The calculation of the contribution of interfering tra- 
jectories of type b significantly more involved. Because 
of the existence of configurations as shown in the bottom 
panel of Fig. [S] in which the encounter region winds one 
or several times around a periodic trajectory, one can- 
not calculate the contribution from trajectories of type 
b by considering two-encounters and three-encounters 
only, as was done in a previous version of this articled 
(Note, however, that since encounters that fully wind 
around a periodic trajectory with period r p exist for pe- 
riods r p < te only, considering two-encounters and three- 
encounters only is sufficient if te < See Ref. 24 for 
a precise calculation that verifies this for closed quantum 
dots.) 

In order to parameterize the combinations of interfer- 
ing trajectories of type b, we use coordinates that mea- 
sure the phase space distance to the periodic trajectory. 
The periodic trajectory, as well as the four interfering 
trajectories of type b are shown again in the top panel of 
Fig- ED The interfering trajectories consist of two "short" 
trajectories (trajectories a.\ and /Ja in Fig. 0] or Fig. 
and two "long" trajectories (trajectories a 2 and (3\ in Fig. 
|3]or Fig. |SJ), where the long trajectories have one more 
revolution around the periodic trajectory than the short 




FIG. 6: Top: Interfering trajectories of type b (dashed and 
solid curves) and the periodic trajectory (dotted). Bottom: 
Short versions of the interfering trajectories together with the 
Poincare surfaces of section. At each Poincare surface of sec- 
tion the stable and unstable phase space coordinates s and u 
parameterize the distance to the periodic trajectory. 



one. We use the two short trajectories as our reference. 
At a point where the phase space distance to the peri- 
odic trajectory is less than the classical cut-off c, for each 
short trajectory i = 1,2 we draw a Poincare surface of 
section and use phase space coordinates (ui, s{), i = 1,2, 
referring to the unstable and stable directions in phase 
space to describe the distance to the periodic trajectory, 
see Fig. Notice that we need to draw two separate 
Poincare surfaces of section because the two short trajec- 
tories do not need to be within a phase space distance c 
from the periodic trajectory at the same time. The long 
versions of the trajectories pass through these Poincare 
surfaces of section twice and have phase space coordi- 
nates (uj, Sje _ATp ) or (uje~ ATp , Si), i= 1,2, where r p is 
the period of the periodic trajectory. 

In order to calculate the action difference AS = S ai + 
S a2 — Sp-L — Sp 2 we perform two successive deforma- 
tions, as shown in Fig. [7| The action difference for the 
first deformation can be calculated using the Poincare 
surface of section drawn in the top left panel of Fig. 
[7| The phase space coordinates of the two trajecto- 
ries are {u2-,s%) and (uie~ A ( Tp ~ r ',sie~ Ar ), where t' is 
the time difference between the two Poincare sections 
in the bottom panel of Fig. El measured along the pe- 
riodic trajectory. Then the corresponding action differ- 
ence is (u 2 -uie- A ( Tp - r '))(s 2 -sie~ AT '), see Refs. l2l29l 
The action difference for the second deformation is cal- 
culated using the Poincare surface of section drawn in 
the right panel of Fig- EI The phase space coordinates of 
the two trajectories involved here are (u\, S2e~ A ( Tp ~ r )) 
and (u2e~ Xr ,si), corresponding to the action difference 
— it2e~ Ar )(s2e~ A<Tp ~ r ) — si). Adding the two action 



8 




FIG. 7: Two successive deformations used to calculate the 
action difference between the pair of trajectories in the top 
left diagram and the pair of trajectories in the bottom left 
diagram. 

differences, one finds 

AS=(u 2 s 2 -u lSl )(l-e- x ^). (22) 

The sign of the action difference is not relevant, as both 
AS and —AS* appear in the final summation. Our ex- 
pression for the action difference AS differs by a factor 
1 — e~ ATp from the action difference used in Ref.|23j This 
difference is unimportant, since relevant periods r p are of 
the order of the mean dwell time td and exp(— Atd) <C 1. 

We count periodic trajectories that consist of several 
revolutions of one shorter trajectory as separate trajec- 
tories. This correctly takes into account the contribu- 
tion from interfering trajectories where the difference be- 
tween interfering trajectories is more than one revolution 
around a periodic trajectory. 

Now we are ready to calculate the contribution B of 
trajectories of type b to the reflection covariance. Re- 
peating the steps used for the calculation of the weak 
localization correction ST, we find 

B NlN " ( FT r^e-^A fdr e-/™ 

l-Tp/2 re e iAS/h-(T s +T u )/T 

x dt 2 I dsiduids 2 du 2 », 9 , . 

J-Tp/2 J-c (^7™) t C nc,l*enc,2 

(23) 

Here t 2 parameterizes the point at which the Poincare 
surface of section for the reflection trajectory of contact 2 
is taken, measured as the time needed to travel between 
the Poincare surfaces of sections for the trajectories 1 
and 2, see Fig. |S| The time i C nc,i is the time during 
which trajectory i remains within a phase space distance 
c from the periodic trajectory, i — 1,2. The division 
by icnc.i cancels a spurious contribution arising from the 
freedom to choose the Poincare surface of section along 




FIG. 8: Definitions of the times appearing in Eq. Xl'-M . The 
top panel shows an example where the two short trajectories 
ai and /?2 (solid curves) have strong correlations before and 
after arriving within the vicinity of the periodic trajectory 
(dotted curve). The time during which these strong classical 
correlations exist is r s or t u . The bottom panel shows an 
example where there are no classical correlations except for 
the time the trajectories are in the vicinity of the periodic 
trajectory. In this case, r s = t u — 0. In the bottom panel, 
the two trajectories have different encounter times t e nc,i and 
t C nc,2- In the example shown in the top panel, the two en- 
counter times tenc,i and t cnc ,2 are (almost) equal. The thin 
solid lines indicate the beginning and the end of the encounter 
with the periodic trajectory; the total encounter regions are 
marked by thick trajectories. The dashed lines indicate the 
positions where the Poincare surfaces of section for the two 
trajectories are taken. 



the trajectory. The times t s and t u indicate the length 
of time over which the trajectories 1 and 2 are correlated 
before and after getting within a phase space distance c 
of the closed loop, respectively, see Fig- El The times Tj, 
j = 1, 2, 3, 4, indicate the duration of the four segments 
of uncorrelated propagation of the two pathes when they 
are outside the encounter region. The classical action 
difference AS is given in Eq. (|22|l above. 
We define the new integration variable 

t' 2 = t 2 - A" 1 ln(c/|ai|) + A" 1 ln(c/|s 2 |), (24) 

which is the time between the points where the trajecto- 
ries 1 and 2 first come within a phase space distance c 
from the periodic trajectory. Without loss of generality, 
we may assume that the phase space coordinate s\ is pos- 
itive. However, we must keep the sign of the phase space 
coordinate s 2 because it enters into the action difference 
AS and into the correlation time t s . Making a variable 
change similar to that of Eq. (JJJ , 

Si = c/<Ti, Ui = cx l a i , i = 1,2, (25) 
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we replace the double integration dsidui, i = 1,2, by 
a single integration dxi, i = 1,2. This variable change 
cancels the denominator i Gn c.iicnc,2 and adds a Jacobian 
(Ac 2 ) 2 , where A is the Lyapunov exponent, see Eq. JSJ. 
We then find 



with 



/•r p /2 ,1 

J( Tp ) = 2A 2 / di 2 / dx x dx 2 J^e-^ 

J-t p /2 J-l ± 

x cos [r(xi =F £2)] i 
where the sign ± is the sign of s 2 and 
c 2 (l - e~ AT p) 



(27) 



(28) 



The prefactor 2 in Eq. I|27[l arises from fixing si to be 
positive. Finally, we change variables U\ = x\, u 2 = ix 2 , 
and w = ±e A ' 2 . The integration over t' 2 and summation 
over the sign ± of s 2 is then represented as the integral 
J dw/X\w\, with e~ Ar p /2 < \w\ < e Ar p /2 , 

f dxiJ f 
J(r p ) = 2X r—r dmdu 2 

J M y_i 

x cos[r(ui - u 2 )]e~ (Ts+T " )/TD . (29) 

The classical limit taken here corresponds to sending r — > 
00, td — * 00, while keeping te/td = (l/ATo)lnr fixed. 
Since / gets multiplied by r 2 ru, see Eq. l|2*5|) above, we 
look for a leading contribution to / of order l/r 2 TD- 

It is instructive to first consider the integral 7(r p ) with- 
out the factor exp(— t„/td). Since t s does not depend on 
the integration variables u\ and u 2 , the integrals over u\ 
and u 2 are straightforward, and one finds 



I(r p ) 



8A . 2 f dw 
— sin r / — < 
H ./ \w\ 



-t b /tt> 



(30) 



This result, however, is deceptive. The fact that sin 2 r — -> 
1/2 for large r is an artifact of our choice of the same 
phase-space cut-off c for the trajectories 1 and 2. The 
same artifact appears when calculating F 1 - 2 for the con- 
tribution of type a: when calculating F\ we discarded a 
rapidly oscillating function oc sinr. However, the square 
of this rapidly oscillating function has a nonzero average 
which should not be retained in the final expression. We 
can avoid this problem by taking slightly different phase- 
space cut-offs for the trajectories 1 and 2, which amounts 
to replacing cos(rwi —ru 2 )) by cos(riUi — r 2 u 2 ). We then 
find 



8A / dw sinn sinr 2 c _ T . /rD 



r\r 2 



(31) 



In the classical limit n, r 2 — > 00 this is a fast oscil- 
lating function and can be discarded. By symmetry, 



replacing the factor exp(— t s /td) by unity also gives a 
fast oscillating function which will be discarded from 
the final answer. Because of this, we may replace each 
factor exp(— t SjU /td) by exp(— t SjU /td) — 1 in the inte- 
grand of Eq. H29|) . After this replacement, we can restore 
r\ = r 2 = r. 



The fact that J(r p ) 



without the factors 



exp(— t s /td) or exp(— t u /td) signals the statistical inde- 
pendence of the conductance and the density of states for 
a chaotic quantum dot with te <C td; — It is also crucial 
to ensure unitarity. Without the factors exp(— t s /td) or 
exp(— t u /td) the expression for I(t p ) contains no refer- 
ence to the phase space distance between the trajectories 

1 and 2 at the point that they approach or leave the pe- 
riodic trajectory, respectively. Without such reference, 
unitarity cannot be preserved, because there will be no 
difference between the case that the trajectories 1 and 

2 originate/terminate at the same contact or not. With 
the factors exp(— t s /td) and exp(— t„/td) unitarity is 
preserved, see the discussion at the end of this section. 

At this point we should specify the times t s and t u . 
The escape of trajectories 1 and 2 is correlated before 
arriving at the closed loop only if the stable phase space 
coordinates si and s 2 have the same sign, i.e., if w > 0. 
In that case, the phase space distance between the two 
trajectories at the point of entrance of trajectory 1 is 
d = c|l — exp(— Xt' 2 )\ = c|l — l/w\ and the phase space 
distance between the two trajectories at the point of en- 
trance of trajectory 2 is d — c\ cxp(Xt' 2 ) — 1| = c\w — 1|. 
Note that for w close to 1 (which is when escape corre- 
lations are relevant) the two phase space distances are 
equal. For definiteness, and in order to have an as- 
signment that is symmetric under the exchange 1 <-> 2, 
we use the largest of these two phase space distances in 
the following considerations. Hence, if w > 1, we have 
d = c(w — 1), whereas d — c(l/w — 1) if < w < 1. 
The correlation time (if positive), then, follows from set- 
ting de XTs = c(b — 1), where e ATp / 2 > b > 1 is a suitably 
chosen number of order unity, 

A-Mn^-lV^-l)] ifl<w<6, 
; A- 1 ln[(&-l)/(l/<t/;-l)] if 1/b < to < 1, 
otherwise. 

(32) 

The final results will be independent of the choice of the 
cut-off b. 

The definition of the correlation time t u is similar. Pro- 
ceeding as before, one finds that t u depends on the prod- 
uct z = wu\/u 2 : this encodes both the time difference be- 
tween the exit points of trajectories 1 and 2 and the sign 
of the unstable phase space coordinates at those points. 
There is one subtlety when regarding escape for the ex- 
iting trajectories: if \z\ < e _ATp / 2 or \z\ > e ATp / 2 , \z\ has 
to be brought back to the range e~ ATp / 2 < \z\ < e ATp / 2 
by multiplication with the appropriate number of factors 
of e ±Ar p. This procedure takes into account that the 
trajectories 1 and 2 can make different numbers of rev- 
olutions around the closed loop before exiting. (There 
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was no such complication when considering r s , because 
the integration domain for w is e~ ATp//2 < \w\ < e ATp / 2 .) 

Since t s and t u appear in the combination 
exp(— T a)U /Tt>) only, we write exp(-T s /r D ) = f(w), 
exp(— t„/td) = f{z). We then have f(x) = 1 if x < 0, 



/(*) = 



a: - 1 
6-1 



1/Atd 



if 1 < x < b, 



and f(x) = 1 if 6 < x < e ATp / 2 . Further, 

f(x) = f(l/x) and f(x) = /(a;e XT »). 



(33) 



(34) 



Note that the function / is continuous and piecewise dif- 
ferentiable. Also notice that / = 1 except in the case 
of strong classical correlations between the trajectories 1 
and 2. 



With this choice of the function f(z), we find that the 
integral I(t p ) reads 



I = — — IfM - 1] / dull dz[f(z) - 1]— sm[ur(z/w - 1)] + / tfa[/(z) - 1]^- sm[ur(zw - 1)] . 
r J w J [J-\w\ az J-i/\w\ az ) 

(35) 

Since f(w) = 1 if w < 1/6 or io > 6, we can restrict the integration range for w to 1/b < w < b. Performing a partial 
integration to z, we have 

I = — — [ b — [/H-l] f^ui /"" dz^MsinM^/™- l)] + / 17 " d z ^Msin[ur(zw-l)]l. (36) 
r 7i/6 ^ Jo \J_ W dz ' y_ 1/l0 <9z 



For the z integration in the second line, there will be contributions from the regions 1/6 < z < w and e~" ATp /6 < 
z < e~ nATp 6, n — 1,2,.... For the z integration in the third line, there will be contributions from the regions 
1/6 < z < 1/w, e~ nXTp /b < z < e~ IlATp 6, n = 1,2, .... In the region 1/6 < z < b one has 

-f(z) if z > 1, 



9z (z — 1)Atd ' 

1 



dz z(z — 1)Atd 
Focusing on the integration range 1/6 < z < 6, we then find 



f(z) ifz<l. (37) 



div f w dz 

I = — — / — / [f(w) — l]f(z) I du {sxa[ru(l / zw — 1)] — sm[ru(z/w — 1)]} 



TD Ji If Ji Z — 1 



(I 



6 du; /""ds, ^ 



+ / / — [f(z)-l]f(w) / du{sin[ru(l/zu;- 1)] + sin [ru(z/u> - 1)1} . (38) 

TTq A «J - 1 Jx Z J 

Although these integrals can be evaluated for general 6, the evaluation is simplest if b — 1 -C 1 (but still 6 — 1 of order 
unity). Writing z — 1 + (6 — 1)C/ U an d tu = 1 + (6 — 1)C/ M an( i expanding in 6 — 1, one finds 

16(6-1) f 1 ^ j d^_ j -[(e/u) 1/ATD - 1]((/m) 1/Atd cos[r(!) - l)e]sin[r(5 - (39) 

I 



rT B Jo JO C Jmax(£,C) U 



Performing the integral over u followed by a partial in- 
tegration over £, the limit r — > oo at fixed ratio te/td = 
(1/Atd) lnr can be taken. One then finds 

/ = 27r 2 (l-r- 2 / ATD )/r 2 T D . (40) 

There is an alternative (and more intuitive) deriva- 



tion of Eq. H39f) if we restrict our attention to interfering 
trajectories that arrive at and depart from the periodic 
trajectory at (classically) close phase space points from 
the very start of the calculation. This is the situation 
drawn in the top panel of Fig. |H1 Instead of choosing 
the Poincare surfaces of section at an arbitrary point 
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during the encounter with the periodic trajectory, we 
may choose them at the beginning and end of the en- 
counter with the periodic trajectory (indicated by the 
thin solid lines in the top panel of Fig. [SJ . At the first 
(entrance) Poincare surface of section, we use the unsta- 
ble phase space coordinate u measured with respect to 
the periodic trajectory, and the difference s' of the sta- 
ble phase space coordinates of the trajectories a± and 
/?2. (Since these trajectories depart from the periodic 
trajectory together, their unstable phase space coordi- 
nate can be considered equal at this Poincare surface of 
section.) Similarly, at the second (exit) Poincare surface 
of section, we use the common stable phase space coor- 
dinate s of both trajectories, measured with respect to 
the periodic trajectory, together with the difference v! 
of the unstable phase space coordinates. The total time 
the trajectory ot\ spends near the periodic trajectory is 
A -1 ln(c/|u|) = A -1 ln(c/|s|), hence \s\ — \u\. In terms of 
these coordinates, exp(— t s /t d ) = |s'/c(6 — 1)| 1 / Atd and 
cxp(— t„/td) = \u'/c{b — l)! 1 /'* 11 " . Fixing the positions 
of the Poincare surfaces of section gives a Jacobian Ac|u| 
and eliminates the factor icnc.i^cnc.2 from the denomina- 
tor in Eq. i|23[l . 2 ™ Upon dividing all coordinates by the 
phase space cut-off c, we then find 



I = A\ 



h-l pi 

ds'du' / duu(\s'/(b - l)! 1 /*™ - l) 

1-b JO 



X (\u'/(b - I^/Atd _ cos [ ur ( s ' _ U /J] j ( 41 ) 

where we added factors 2 for the signs of u and s. Rewrit- 
ing the integrals such that the integrations over s' and 
v! are all between and 6—1 and performing a partial 
integration to s' we find 



6-1 



ds' 



T du' / du(\u'/(b- 1)\V Xt ° - 1) 



J 

x \s'/(b - l)) 1 /^ s i n ( ws ') cos(wru'). (42) 

Upon shifting variables s' = £(6—l)/u and it' = £(6— l)/u 
one then arrives at Eq. above. 

The z integrations with 1/6 < ze nrX < b with n = 1, 
2, ... do not give a contribution to I(t p ) in the limit 
r — > oo. This can be seen by noting that all oscillating 
integrals all contain fast oscillating phases proportional 



to r(l 



-n\r-r 



Putting everything together, we find 



B = 



N?N§ 



2{N 1 +N 2 ) i T D J p (l 



rfr n 



-t p /t 



(1 



3 -2t e /t d \ 



(43) 

Setting a lower cut-off for the r p -integration at r p > 1/ A 
and taking the limit tdA 3> 1, we finally arrive at the 
simple result 



B = 



"(1 



-2t e /td 



2(JVi + N 2 y 

Substitution into Eq. I|20|) then gives 

„2 \ 2 *r2 AT2 



varf; = 2 ( ^~ 
h 



(N r + N 2 ) 4 



(44) 



(45) 



Equation (|45|l is the main result of this article. The 
variance of the conductance is independent of the Ehren- 
fest time. Equation lj45|l was derived for a quantum 
dot without an applied magnetic field. With a magnetic 
field strong enough to fully break time-reversal symmetry 
varG is reduced by a factor two. 

The most remarkable feature of Eq. (|45|l is that the 
conductance fluctuations survive in the limit te 3> td- 
The particular classical trajectories that give rise to con- 
ductance fluctuations in this limit can be identified by 
inspection of the calculation above. The constant term 
in Eq. I|45|l can be traced back to the lower limit on the 
u integration in Eq. I|39f) which, in turn, results from tra- 
jectories that wind many times around the periodic tra- 
jectory. Although such trajectories must spend a time 
> te inside the quantum dot in order to contribute to 
the conductance fluctuations, their survival probability 
depends on the period r p of the periodic trajectory only, 
not on the actual time they spend inside the quantum 
dot. 

Instead of calculating the conductance variance 
through the covariance of the reflection off the two con- 
tacts, one can also directly calculate the variance of the 
conductance or the variance of the reflection. As in the 
case of the calculation of the quantum correction to the 
average conductance, there are two types of contributions 
to the fluctuations: interfering trajectories for which all 
encounters lie within the interior of the quantum dot and 
interfering trajectories for which at least one encounter 
touches the lead opening. The calculation of the first 
type, all encounters inside the quantum dot, proceeds in 
essentially the same was as the calculation of the reflec- 
tion covariance outlined above. Below, we discuss how 
the above calculations should be modified to include the 
second type, for which one or more encounters touch 
the lead opening. We find that, once encounters that 
touch the lead opening are taken into account, unitarity 
is obeyed for contributions of type a and b separately. 

For trajectories of type a, there are two encounters that 
each can be close to a lead opening. For each encounter, 
however, the configuration of trajectories is identical to 
that of coherent backscattering calculation. Repeating 
the steps of the last part of the previous section, verifi- 
cation of unitarity is immediate. 

For trajectories of type b, a schematic drawing showing 
the difference between an encounter that fully resides in- 
side the quantum dot and an encounter that touches the 
lead opening is shown in Fig. For an encounter that 
touches the lead opening, one replaces the phase space 
coordinate differences s' or u' in Eq. I|41|l by t>, where v 
represents the perpendicular component of the momen- 
tum in the lead opening. As in Sec. ^ we normalize 
v such that the volume element in phase space is dudv. 
With this choice of phase space coordinates, the action 
difference for two pairs of trajectories involved in an en- 
counter that resides inside the quantum dot is the same 
as the action difference for two pairs of trajectories for 
which one of the encounter extends to the lead opening, 
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FIG. 9: Comparison of a part of a type b encounter that fully 
resides inside the quantum dot (left) and an encounter that 
touches to lead opening (right). The relevant segments of the 
(short) trajectories ai and P2 are shown solid, the periodic 
trajectory is shown dotted. 



ference between weak localization and conductance fluc- 
tuations not only concerned the dependence of their mag- 
nitude on the Ehrenfest time — exponential decay versus 
independence of the Ehrenfest time — , but also the mini- 
mal dwell time of trajectories that contribute to quantum 
interference. According to the numerical simulations, the 
minimal dwell time for weak localization is 2te, whereas 
the minimal dwell time of trajectories contributing to 
conductance fluctuations was found to be te- Below we 
show that this observation is fully consistent with the 
semiclassical theory. 

We first consider weak localization. With the imagi- 
nary term added to the energy, one finds that the quan- 
tum correction to the transmission is given by 



up to the substitution s' — > v or u' — » v. One then ob- 
tains the contribution of encounters that touch the lead 
opening by making the replacement 



Ni 



Nx+N 2 \ 1L I r D 

\.7=1 



n . ^ 



,-Tj/td 



,/ru 



T D 



(46) 



in the expression for the correlator with encounters that 
reside inside the quantum dot only. Here r' is the travel 
time between the lead opening and the point where the 
two trajectories get close to the periodic trajectory and 
exp(— Tu/td) = v/c(fo — 1) | 1 / A " r D _ The r' integration gives 
a factor 1 — e _Tl, / TD , so that the remaining integrals are 
the same as for the case that all encounters are inside the 
quantum dot and one verifies that unitarity is obeyed for 
type b encounters as well. 



IV. TIME DEPENDENCE 

Although the calculation of the previous section an- 
swers the question which trajectories are responsible for 
the conductance fluctuations in the limit of large Ehren- 
fest times — trajectories that wind many times around a 
certain periodic trajectory — , it does not tell us how long 
these trajectories spend inside the quantum dot. That 
question can be answered by adding an imaginary term 
to the energy. Such an imaginary term gives rise to an 
additional exponential decay exp(— i a /2T a b s ) for each tra- 
jectory a, where t a is the total duration of trajectory a 
and T a bs is the corresponding absorption time. A minimal 
time needed for quantum interference shows up through 
an exponential dependence on l/r a b s . In the calculations 
below, we keep the ratios te/td and t-e/t^s fixed while 
taking the classical limit h — » 0. 

The addition of absorptio n has b een used in the nu- 
merical simulations of Refs. j_6|3C|.34| as a diagnostic tool 
to investigate the microscopic mechanism of quantum in- 
terference corrections. Reference 30 found that the dif- 



5T = 



AxA 2 



dr., 



n . % 



dT,e- T3/TD - T3/T 



x / dsdu- 



isu/ft — tonc/Tt>— 2t em! /-7V 



2nht c 



(47) 



Note the factor two in front of £ C nc/ r abs m the survival 
probability during the encounter. This factor two arises 
because trajectories that contribute to weak localization 
travel the encounter region twice. Performing the inte- 
grations as in Sec. |nl one finds 



ST = (1 + 2^3/7-^) ^ 

— TE/TD— 2TE/T al 



(Wi+Wa + Wrabs) 3 
Similarly, for reflection one finds 



(48) 



SRt 



N,N 2 (l + 2r D /r abs ) c _ rE/TD _ 2TE/rabs 
(N 1+ N 2 f (l + r D /r abs ) 3 

. Ni (W^abs) 2 „-r E /r n -2r E /r^ 

Ni+N 2 (i + TD / Tahs f 



(49) 

In the limit te <td, r a b s these results agree with those 
obtained from random matrix theory^ Note that SRi no 
longer equals — 5T, i — 1, 2, in the presence of absorption. 

The exponential dependence oc exp(— 2rE/r a b s ) indi- 
cates that the minimal duration of a trajectory that con- 
tributes to weak localization is 2te: a self-encounter lasts 
one Ehrenfest time and each trajectory to weak localiza- 
tion passes through the encounter region twice. The fact 
that the dwell time dependence oc exp(— te/td) does not 
have this factor two is a result of classical correlations 
between the two segments of the trajectory that pass 
through the encounter regioni&S, Both effects — the min- 
imal time 2te needed for weak localization and the expo- 
nential suppression oc exp ( — te / t n) are consistent with 
the numerical simulations of Refs. 1(11.3(1 

For the conductance fluctuations we treat the contri- 
butions of trajectories of types a and b separately. We 
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limit our discussion to the covariance cov(i?i, R 2 ) of the 
reflections off the two point contacts. Including absorp- 
tion into the covariance contribution A of trajectories of 
type a, one finds 



A 



= (A^ 2 ) 2 (1 + 2r D /r abs ) 2 TE/TD _ 4TE/rnha 
(iVi+iV 2 ) 4 (l + r D /r abs )6 



(50) 

In order to calculate the contribution of trajectories 
of type b, we need to calculate the integral I(t p ) in the 
presence of absorption, 

J(t p ) - 2\e- T ?^ fp. [ d Ul du 2 \u lU2 \ 1/XT ^ 
J \ w \ J-i 

x cos(ntti - r 2 u 2 )e- (Ts+T " )(1/rD+2/Ti > bs) . 

(51) 

Here we used that the time the short trajectories spend 
near the periodic trajectory is A -1 ln(l/|ui|), i = 1,2, 
and that the time spent by the long trajectory is r p + 
A -1 ln(l/|ztj|). As before, we first calculate I without 
the exponential factors involving r s or t u . We then find 



I(r P ) 



Att 2 t, 



£.g— T p /r a bB— 2te /r a bB 



(52) 



Replacing one exponential factor exp[— t SiU (1/tu + 
2/r abs )] by exp[-r s , u (l/r D + 2/r abs )] - 1 while still set- 
ting t u>s — in the other exponential factor, one finds no 
significant contribution to /. Hence, the remaining con- 
tribution to I can be calculated by replacing both factors 
exp[-T s ,„(l/TD+2/T abs )] by exp[— r s , u (l/rD+2/r a b s )]— 1. 
The result then follows from Eq. l|39ll after addition of a 
factor u 2 / XTabs and replacement of 1/tq by l/rD + 2/T abs . 
Performing the integrals as described in Sec. IIIII we find 



'(Tp) 



9 2 * 



T v /r a ba -2r E /r a ba 



27r 2 (l/T D + 2/r abs ) ; 
r 2 (l/TD + l/r abs ) 



x e 



D -2TE/T aba _ „-2r E /r D -4T E /T ab s 



(53) 



Adding the contributions of trajectories of types a and 
b, the final result becomes 



cov(i?i, R 2 ) 



(JV1JV2) 2 



= -2r E /Taba 



(54) 



{N!+N 2 y 

(1 + 2r D /T abs ) 2 + 2(r D /r abs ) 2 
(l+T D /r abs ) 6 

The limit te/td — ► agrees with the result obtained 
from random matrix theory™ The exponential decay 
tx exp(— 2te/td) corresponds to a minimal dwell time 
te for trajectories that contribute to conductance fluctu- 
ations. This is consistent with the numerical simulations 
of Refs.|jy] and|13- The same conclusion remains true for 
other correlators of transmissions and reflections (which 
need not be equal to cov (R±, R 2 ) in the presence of ab- 
sorption). 



V. CONCLUSION 

In the previous sections, we have calculated the 
Ehrcnfest-time dependence of the weak localization cor- 
rection and the conductance fluctuations of a ballis- 
tic quantum dot with ideal point contacts, using a 
trajectory-based semiclassical theory. We find that the 
Ehrcnfest-time dependences of weak localization and con- 
ductance fluctuations are remarkably different: whereas 
weak localization is suppressed exponentially if the 
Ehrenfest time te is much larger than the mean dwell 
time rrj, the conductance fluctuations are independent 
of the ratio te/tq. Our calculation explains the numer- 
ical simulations of the conductance fluctuations in Refs. 

mm 

The numerical observation of Ehrenfest-time indepen- 
dent conductance fluctuations was remarkable, because 
trajectories need to remain inside the quantum dot dur- 
ing at least a time te if they are to contribute to quantum 
interference corrections, so that one expects that inter- 
ference corrections to the conductance should disappear 
if te ^> td • Our calculation shows that the latter expec- 
tation is not always born out. The Ehrenfest-time inde- 
pendent contribution to the conductance fluctuations we 
find here arises from pairs of trajectories that wind many 
times around a periodic trajectory. On the one hand, 
such configurations of interfering trajectories have suffi- 
cient duration to allow for small action differences and 
thus provide a significant quantum interference correc- 
tion. On the other hand, they are not as much affected 
by classical escape into the leads: The survival proba- 
bility depends on the period of the periodic trajectory 
involved, not on the Ehrenfest time. 

In Refs. ^3 an d E] the 'effective random matrix the- 
ory' of Silvestrov et aL 36 was used to explain the nu- 
merical observation of Ehrenfest-time independent con- 
ductance fluctuations. According to the effective ran- 
dom matrix theory, phase space is separated into a 'clas- 
sical part', corresponding to all trajectories with dwell 
time less than the Ehrenfest time, and a 'quantum part', 
which has all trajectories with dwell time larger than 
te- The wave nature of the electrons plays no role in 
the classical part of phase space, whereas the quantum 
part of phase space is described using random matrix 
theory. The effective random matrix theory was able to 
correctly describe the Ehrenfest-time dependence of shot 
noised and the density of states of an Andreev quantum 
dotj2£i2£iSi but it missed the Ehrenfest-time dependence 
of weak localization^ According to our semiclassical the- 
ory, the effective random matrix theory not only gave the 
correct magnitude of the universal conductance fluctua- 
tions, but it also gives the right dependence on an imagi- 
nary potential. 34 While it is understood that the effective 
random matrix theory is not a comprehensive theory of 
the Ehrenfest time dependence of quantum transport, it 
remains an interesting question to classify which trans- 
port phenomena are described by this phenomenological 
theory and which are not. 
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